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ABSTRACT 

In this paper we present a new method for selecting blue horizontal branch 
(BHB) candidates based on color-color photometry. We make use of the Sloan 
Digital Sky Survey z band as a surface gravity indicator and show its value for 
selecting BHB stars from quasars, white dwarfs and main sequence A type stars. 
Using the g, r, i, and z bands, we demonstrate that extraction accuracies on par 
with more traditional u, g, and r photometric selection methods may be achieved. 
We also show that the completeness necessary to probe major Galactic structure 
may be maintained. Our new method allows us to efficiently select BHB stars 
from photometric sky surveys that do not include a u band filter such as the 
Panoramic Survey Telescope and Rapid Response System. 

Subject headings: Stars: horizontal-branch - Techniques: photometric - Techniques: 
spectroscopic 
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1. Introduction 

Blue horizontal branch stars (BHBs) are helium core-burning stars, typically about 
0.7Mm and 3R^, w hich exist at a near-constant absolute g band magnitude of about 0.7 



( lYanny et al 



20001 ). Since they are such luminous objects with a predictable brightness, 
they are frequently used as standard candles to explore distant Galactic structure and 
are widely sought tracers for structural and kinematic studies of the Milky Way (e.g. 
Clewley & Jarvis 2006, Ruhland 2011, Wilhelm et al. 1999, Xue et al. 2008). Owing to 
their old ages (>12 Gyr in many globular clusters, Dotter et al. 2010), BHB stars are 
particularly excellent halo tracers, since the halo is expected to be an old, metal-poor 
structure. Other high luminosity stars, such as K and M giants, are less desirable as their 
absolute magnitudes cover a wide range of values and they may intrinsically be much 
younger objects. 

Photometrically, BHBs reside in a specific color range, being just bluer in g — r than 
the instability strip (RR Lyraes) and just redder than extended vertical blue horizontal 
branch (sub-dwarf) stars. This "blue" color slice (-0.3 < g — r < 0.0; all magnitudes and 
colors in this paper are extinction corrected and dereddened according to the values of 
Schlegel et al. 1998 unless otherwise stated) contains three main contaminants to BHB 
selection: distant quasars, foreground white dwarfs, and main sequence A stars (MSA; blue 
stragglers are photometrically very similar to MSA stars and so they are included in the 
same category for the purposes of this paper). 

The most accurate way to separate a BHB star from these contaminants is via 
spectroscopy. Spectral templates will immediately eliminate white dwarf and quasar 
contamination (as white dwarfs have very broad absorption lines, and quasars have strong 
emission rather than absorption lines). To separate MSA stars from BHBs, two main 
approaches exist: first a comparison of Balmer line depth (f m ) and broadening (-D0.2) 
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caused by the different surface temperatures and surface gravities, respectiv ely (see Pier 



'. 983) , of the stars; and secondly the scale-width-shape method described by 



Clewley et al 



( 120021 ) which separates the two species based on the Balmer line fits to a Sersic profile. 



In the absence of spectroscopic data however, photometric methods may be applied to 
separate BHB stars from othe r blue species . This has been shown to be quite effective in 



several studies. For example: 



Yanny et al 



( 120001 ) demonstrated a filter cut in the Sloan 
Digital Sky Survey (SDSS; York et al. 2000) which separates BHB stars from MSA stars 
adequately enoug h to discern significant structure in the plane of t he celestial equator ; 



Sirko et al 



Yanny et al 



(120001 ) on 



(|2004[ ) use a "stringent" cut similar to that proposed by 
the basis of spectroscopic data and a combination of the scale-width-shape method and 
the Dn.<? and f m methods in their kinematic studies of the Galactic halo; more recently 



Bell et al. 



( 120101 ) further refine this "stringent" cut in their investigation of the ratios of 



BHB to main sequence turn off stars in the halo. 



For the SDSS, 



Lenz et al. 



( 119981 ). using Kurucz model spectra, suggest that optimal 



gravitational separation for blue A-type stars lies in the u — g color space. This is largely 



due to t 
gravity. 



re Balmer jump, a spectral feature at 365 nm which is highly dependent on surface 



Lenz et al. 



( 119981 ) also suggest separation in i — z space, but note that gravitational 
splitting in this color-space is less effective than in the u — g color space. We show that 
splitting here is caused by the Paschen features which reside in the z band and are also 
sensitive to surface gravity. 

This gravitational splitting in the z band is a serendipitous development with regards 
to the ongoing Panoramic Survey Telescope and Rapid Response System (Pan-STARRS; 
Kaiser et al. 2010) project, which features filters similar to the SDSS photometric system 
(ugriz, see Fukugita et al. 1996, Gunn et al. 1998). The main difference being that the 
SDSS, originally designed as a primarily extragalactic survey, uses the u band to efficiently 
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separate out quasars (e.g. Richards et al 2009 and references therein) while Pan-STARRS 
is primarily concerned with near earth objects and, as such, opts for a near infrared y 
band instead of the u filter. It is worth noting that Pan-STARRS, a purely photometric 
survey, plans to su rvey three-fourths of the celestial sphere to 24 th magnitude in the optical 



(iKaiser et al. 



20101 ) while SDSS, as of Data Release 8 (DR8; Aihara et al. 2011) has imaged 
one-third of the sky to a limiting magnitude of 22.2 in g and r. So, by developing an efficient 
way to select BHB candidate stars in the absence of a u band, we may facilitate statistical 
structural studies in a larger area of sky, to fainter magnitudes, than is currently possible 
with SDSS data. 

With this in mind, we create a color space cut to select out BHB stars from the 
primary contaminants based on surface gravity measurements, but rather than exploiting 
u band separation, we explore the usefulness of the z band. With an end-goal of structure 
mapping and halo studies in mind, we aim primarily for a high purity selectio n at the cost 



of com plet eness. We attemp t to achieve accuracies similar to those quoted by 



Sirko et al. 



(12004] ) and 



Bell et al.l (120101 ) (< 30% contamination), however, we also expect a lowered 



sample completeness because of the poorer gravitational separation in this color space. In 
§2 we describe the formulation of the color cut using SDSS spectroscopic data. In §3 we 
estimate the completeness of the cut by examining a sample of 10 globular clusters located 
in the SDSS footprint. In §4 we explore the usefulness of this method for probing Galactic 
structure by examining photometric data of the celestial equator which is rife with features. 
In §5 we discuss the implications of this method. 



2. The Color Cut 



To investigate t he separation of various blue objects in the i — z color space suggested 
by iLenz et al.l (119981 ). we select spectra from SDSS DR8. DR8 extends the SDSS footprint 



-6 - 



to now cover a full third of the celestial sphere and increases the total number of spectra to 
over 1.8 milHorJ^. 

We selected spectroscopic data from the entire SDSS footprint, only clean (not near 
saturated pixels) objects evaluated as point sources were selected. We used only primary 
measurements (in the case of multiply observed objects, the highest signal to noise readi ng 



i s flag ged as the primary on e). To make our study consistent with those of 



(2004 


) and 


Bell et al. 


(2010) 



Sirko et al. 



( 120101 ) we select only stars with g < 18. 



In Figure 1, we plot the adopted log(g) value as calculated by the SEGUE Stellar 
Parameters Pipeline (SSPP, see Lee et al. 2008; SEGUE stands for the "Sloan Extension 
for Galactic Understanding and Exploration", see Yanny et al. 2009) against g — r. We 
note that in the blue color space of -0.3 < g — r < 0.0 there is excellent differentiation of 
the lower gravity BHB stars, residing in the range of 3.0 < log(g) < 3.75, and the higher 
gravity MSA stars (3.75 < log(g) < 5.0). These color and gravity ranges form the basis for 
our identification of these two types of star. 

Some contaminants have no SSPP atmospheric parameters, so we cannot identify 
them based on their log(g) estimates: the two principal examples are white dwarf stars 
and quasars. White dwarf stars raise critical flags in the SSPP because of the width of 
their Balmer li nes (Do ■? > 35. A) and quasars raise critical flags due to their strong 



emission lines ( ILee et al 



20081 ). To identify these contaminants, we rely on the ELODIE 
template matches (the ELODIE archive is a set of high resolution spectroscopic readings 
collected using the ELODIE spectrograph, which has been operating on the Observatoire 
de Haute-Provence 1.93 m telescope since 1993; Moultaka et al. 2004) as output by the 
SSPP. The remaining contaminants are binned together- they mostly consist of A and F 



1 http:/ /www. sdss3.org/dr8/ 
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stars for which the SSPP gravity reading was either inconclusive or outside the bounds of 
the prior mentioned BHB and MSA star bins. 

Using this data set we look to construct a photometric color cut in g — r (a temperature 
indicator) vs i — z (a surface gravity indicator) color space (see Figure 2). To do this we use 
a k- nearest-neighbors classification algorithm. This algorithm classifies unknown objects 
based on their proximity to known objects. We compare a uniform grid in g — r vs i — z 
color space (the grid consisting of 101x101 nodes over the color space -0.3 < g — r < 0.0 
and -0.25 < % — z < 0.05) to the spectroscopic data- grid nodes with 50% of their 5 nearest 
neighbors having been classified as BHB stars were said to reside in BHB color space. 
We then drew a rough selection box around this BHB color space defined by the points: 
(g-r,i-z) = (-0.30,-0.18), (-0.05,-0.06), (-0.05,-0.02), (-0.30,-0.02). 

In Figure 2 we plot these 5 datasets: high gravity MSA stars, low gravity BHB 
stars, spectroscopically identified white dwarfs, spectroscopically identified quasars and 
inconclusive points. The lower panel of the plot is a g — r vs i — z color-color plot: BHB 
stars are plotted with diamonds while contaminants are plotted with points. Note that 
only 1 in 5 data-points are shown in the lower panel of the plot to avoid obscuration. 
However, the upper histogram shows all of the objects which are inside the selection box. 
It is apparent that this separation culls white dwarfs with acceptable efficiency. Note that 
type DC white dwarfs were neglected from this analysis due to a lack of a SSPP "DC 
white dwarf classification and their intrinsically low numbers (<4% of all white dwarfs). 
Additionally, one could use a proper motion diagram to remove some of the remaining 
white dwarfs and some of the closer MSA stars. 

We do not get much more separation from the stellar contaminants in other color 
spaces. However, we can more efficiently select out quasar contaminants in the g — r vs 
g — z color space. We use g — z color space due to the spectral profiles of blue stars which 
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are characterized by the tail end of blackbody profiles (tending to lower g — z values) and 
quasars which are more uniform emitters (tending to higher g — z values- especially in the 
case of high redshifts). In Figure 3 we present this additional color cut. This cut is defined 
by the points: (g-r,g-z) = (-0.3,-0.72) , (-0.3,-0.57) , (-0.05,-0.08) , (-0.05,-0.23). We plot 
only the data which pass the cut shown in Figure 2. Quasars are plotted as open circles 
and stellar sources are plotted as dots to accentuate the separation. Again histogrammed 
above are the objects passing the color cut. As expected, we see very little change in the 
stellar contaminants, but quasar contamination drops significantly. 

As a test of accuracy, we consider the number of stars falling within the color cut 
which are not BHB stars: out of ~4300 spectra which pass this cut, we find ~77% 
of them to be BHB stars; this represents a sample of ~74% of the stars originally 
identified as BHB stars. F or comparison, w e run a similar test on the same data using 



the color cut suggested 



employed by 



^y 



Sirko et al. 



Bell et al 



(120 lOf ) which is similar to the "stringent" color cut 
f l2004h (0.98 <u- g < 1.28 , -0.2 < g-r < -0.06 , excluding 
{[u-g- 0.98]/0.215) 2 + {[g - r + 0.06]/0.17) 2 < 1). Using this refined u band dependent 
test, we similarly select a sample that is ~74% pure and ~72% complete. This test is biased 
by the SDSS spectroscopic selection algorithms, and so is neither a strict test of purity nor 



completeness. 



3. The Globular Cluster Test 



We use a comple teness test similar to the one employed by iBell et al. 



(120 loh . Using 



Jordi fc Grebell ( 120101 ) as a reference, we select 10 globular clusters from the SDSS footprint 
with pronounced BHBs. By running the constituent stars through the photometric cut 
described above, we may get a second measure of how effective the algorithm is at extracting 



BHB stars. The clusters chosen and pertinent information is tabulated in Table 1. 
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Since the SDSS photoObj All pipeline fails for crowded fields, the cores of dense 
objects, such as globular clusters, are often omitted from the general photometric data. For 
accurate and complete cluster photometry we turn to the SDSS "value added" catalogj^- 
in particular we use the "ugriz DA OPHOT photom etry of SDSS+SEGUE Globular and 



Open Clusters" catalog produced by 
Since the 



Anetal 



(}2008|). 



An et al 



(120081 ) value added catalog was created with SDSS Data Release 
Seven (DR7; Abazajian et al. 2009) imaging and pipelines, we use DR7 photometric 
data for this cluster analysis. Data are selected from DR7 within one third of the tidal 
radii of the relevant clusters (Harris (1996) Catalog, 2010 revision) and then matched to 
the DAOPHOT photometry. In the cases of duplicate points, DAOPHOT photometry 
is preferred. The stellar magnitudes and colors are extinction corrected and dereddened 
according to the Harris Catalog E(B-V) values. 

By examining the color-magnitude diagrams of the clusters, we visually select 
1-magnitude wide boxes encompassing the main portion of the blue horizontal branch in 
the -0.3 < g — r < 0.0 regime. This selection box contains the stars we will consider "true" 
BHBs. All stars outside of these boxes will be considered "false." It is apparent from 
clusters such as NGC 7078 and NGC 5272 that these 1-magnitude boxes suffer RR Lyrae 
contamination on the red end and in the clusters NGC 6205 and NGC 6341 we see hot 
sub-dwarf contamination on the blue end. However, in general, these boxes should consist 
mainly of BHB stars. The ranges for these magnitude boxes are also given in Table 1 and 
the cluster color-magnitude diagrams with the selection boxes superimposed are shown in 
Figure 4. 

We select only stars passing the color-cut selection described in §2. We do not consider 



2 http:/ /www. sdss.org/DR7/products/value_added/index. html 
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this test to be a test of accuracy due to the selective enrichment of BHB stars in these 
cluster-fields- however, we do consider this to be a good test of completeness. Stars passing 
our z based color cut constitute a sample that is ~95% pure and ~51% complete. Similar 
to §2, we cut the data using the u based "stringent" color selection and find a sample that 
is ~92% pure and ~57% complete (see Table 2). 

We consider the statistics on the spectroscopic data to be a better indicator of accuracy 
and the globular clusters to be a better indicator of completeness. Thus, when considering 
this test in conjunction with the spectroscopic statistics, we consider the u based cut to be 
74% pure and ~57% c omp 



(120 101 ) and 



Sirko et al. 



ete which is in agreement with the values quoted by 



Bell et al. 



(120041 ) . The z based cut similarly has a ~77% purity and a ~51% 



completeness. 



4. The Celestial Equator 



To show the ability of this selection method to probe features of the Milky Way's 
structure, we examine the celestial equator. This area (-1.26° < declination < 1.26°) has 
been almost completely imaged in the SDSS, excepting areas too close to the Galactic 
plane. The portion of t his great circle above the Galactic equator is known as stripe 10 (see 



Stoughton et al 



(120021 ) for stripe naming conventions in the SDSS) and the portion residing 
in the southern Galactic hemisphere is stripe 82. Stripe 82 has been imaged repeatedly to 



20081 ). RR Lyraes 



promote studies of variable objects, such as supernovae (IFrieman et al. 
and quasars. As such, this section of the sky has be e n ext ensively probed for evidence of 



structure (see 



Newberg et al 



(120021 ) and 



Sesaretal 



(120071 ) for examples) 



We create our sample by selecting all photometrically clean sources identified as point 
sources (selecting only primary measurements in the case of duplicates) residing in the 
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plane of the celestial equator and pass their dereddened photometric data through the 
color cut suggested in §2. We note that, in general, noise in spectra of stars fainter than 
g—18 throw errors into the parameter determinations (Sirko et al. 2004; exceptions include 
SEGUE pencil beams which were imaged for various time scales to maximize signal to noise 
ratios at two main magnitude bins, see Yanny et al. 2009) - but since the CC D camera can 



19981 ) . photometric 



reliably determine every color to much fainter magnitudes (IGunn et al. 
separation is still practical. Photometric separation has the parallel benefit of being viable 
for a much more complete sample of stars when compared to spectroscopic methods. 

The results of our photometric selection are plotted in Figures 5 (the northern Galactic 
hemisphere portion of the celestial equator) and 6 (the southern Galactic hemisphere 
along the celestial equator). Immediately evident are several well documented structures: 
among these are the Sagittariu s stream in both the no rthern and southe rn Galactic 



hemispheres (I Yanny et al. 



20001 ). the Virgo overdensity (IVivas et al 



20011 ) in the north, 



and the Hercul e s-Aqu ila cl oud in both the north and the south noted independently by 



Newberg et al. 


f 


2002 


) and 


Belokurov et al. 


(2007 


), and the anomalous density at R.A. = 


160° noted by 


Newbere et al. 


(2002 


)• 



5. Spectral Interpretation 



We suspect that the z band shows gravitational separation due to the gravity sensitive 
Paschen features (the Paschen analog to the Balmer jump resides at 8201 A), which lie in the 
z band. We anticipate that this separation may be even more effective on the Pan-STARRS 



telescope which has a much more in 
and thus a more responsive z filter 



rared sensitive CC D than the SDSS imaging camera, 



Stubbs et al. 



2010). Figure 7 depicts spectra of the 



objects discussed in this paper (BHB stars, MSA stars, white dwarfs, and quasars) and 
includes the throughput functions of both the SDSS 2.5 m telescope and the Pan-STARRS 
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1.8 m telescope. 

To explain the difference in % — z colors of BHB stars and MSA stars, we examine 
high signal-to-noise spectra from the SDSS survey with particular attention paid to the 
7000A-9200A range. We collect FITS images of all spectra in the SDSS survey which have 
a signal-to-noise ratio > 50, are primary science readings, and classified as stars. Similar to 
§2, we select stars in the -0.3 < g — r < 0.0 color range and define BHB stars as those in 
the 3.0 < log(g) < 3.75 gravity range; MSA stars are defined as having 3.75 < log(g) < 5.0. 

We construct two "super-spectra" by combining all spectra from each population. The 
spectra were individually shifted in wavelength such that the minimum of their Balmer-a 
absorption feature fell at 6563A to account for differing radial velocities. They were then 
normalized to their flux value at 7500 A. It became apparent that, since our sample consisted 
of the whole -0.3 < g — r < 0.0 color range (which spans from about 7500K to 10000K) 
a systematic difference in temperature needed to be accounted for (as our BHB sample is 
centered on this color range, but our MSA sample trends toward cooler, redder colors). 
We select BHB stars in the color range -0.15 < g — r < -0.10 and MSA stars in the color 
range -0.17 < g — r < -0.12- these color ranges produced composite spectra which had 
approximately equal continua and Balmer line depths. We then bin the spectra into 2A 
wide bins and accept the median values of these bins as the composite value. 

Figure 8 shows the well known Balmer-a absorption feature for our two spectra. We 
easily discern the differing profile widths and shapes that are the basis for differentiation 
of these types of stars in spectroscopic studies (e.g. Pier 1983, Clewley 2002; the BHB 
features are significantly slimmer than their MSA analogs). In the second frame we expand 
the Paschen features of these two types of stars. It is unsurprising that we see a similar 
effect in the z band features- the BHB features are significantly slimmer. Near a hydrogen 
absorption series limit, the features will overlap to form a quasi-continuum, thus the more 
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intense Stark pressure broadening of the MSA stars makes their features wider and forces 
their pseudo-continua lower. This combination of slimmer absorption features and a higher 
Paschen continuum explains the difference in z band magnitude for BHB and MSA stars at 
the same temperature. 

An interesting corrollary to this result is the possibility of using only the Paschen 
features to separate BHB and MSA stars spectroscopically. To examine this, we consider 
only the area of the spectra between 8500Aand 9000A; we only use this portion of the 
spectra since fitting to the entire spectra provides especially poor fits in the Paschen region. 
The continuum is defined to be the portion of the spectra exactly in between any two 
Paschen minima ±5A (2 standard deviation outliers are discarded)- the continuum is 



then fit to a terti atic, which is divided out. The resulting normalized spectra are analyzed 



using an IRA1 



(llbdyj 



19861 ) fitting routine. Included in Figure 8 is a plot of intensity vs 
the Voigt profile equivalent width of the Paschen line at ~8596A (the Paschen 14 line) for 
the entire high signal to noise dataset. We see good separation of BHB and MSA stars in 
this plot. Notably, in this region, we see the BHB stars as having larger equivalent widths, 



which is counterintuitive. This is an effect of the pseudo-con tinuum shrin 



features artificially owing to their broad and extended wings ( iFremat et ai- 



ring t he MSA 



3) 



1996)- an effect 



corroborated by the correspondingly higher core intensities of the BHB stars. 



3 IRAF is distributed by the National Optical Astronomy Observatories, which are oper- 
ated by the Association of Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 
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Conclusion 



In §3 we find, in agreement with iBell et all (120101 ) and lSirko et all ( 120041 ). that u based 



color cuts may photometrically select out samples of BHB stars which are ~74% pure and 
~57% complete. We also show that the unexplored z band is capable of selecting out 
samples of BHB stars with similar accuracies (for structural studies, purity is more valuable 
than completeness), having achieved an estimated purity of ~77% based on spectroscopic 
data and globular cluster analysis. However, since the A-star separation by gravity is 
weaker in the z band than in the u band (as noted by Lenz et al. 1998), we obtained a 
lower completeness of ~51%. For increased accuracy, we suggest a combination of the two 
techniques, since they extract overlapping but non-identical sets. 

In §4 we show that this color cut allows structure mapping in the Milky Way to large 
distances. Structures such as the northern Galactic portion of the Sagittarius tidal stream, 
which are on the very edge of F-turnoff detectability, are easily discerned in their radial 
completeness. Thus this method may be used for large-scale structural analysis in addition 
to high efficiency spectroscopic targeting. In §5 we show the spectroscopic reasoning for 
the z band magnitude differences in BHB and MSA stars- namely that Stark pressure 
broadening causes the continua of these stars to exist at fundamentally different levels. 
In short, we suggest that this algorithm is an excellent addition to existing photometric 
selection methods for structure characterization via BHB stars. 
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Fig. 1. — Dereddened g — r color versus log(g) as calculated by the SSPP pipeline for all 
spectroscopic sources with valid log(g) measurements. In the A-star colors (ranging from 
about —0.3 < g — r < 0.0) we see two distinct stellar species: the lower surface gravity blue 
horizontal branch stars and the high surface gravity MSA stars. It is important to note that 
white dwarfs are largely absent in thi s diagram (as ar e quasars) since the SSPP adopted 



log(g) algorithm fails for these objects. 



Leeetal 



(120081 ) estimate the average uncertainty in 



the log(g) measurement from the SSPP to be 0.29 dex. 
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Fig. 2. — The lower panel of this plot shows the distribution of objects in the g — r versus 
% — z color space. Plotted with diamonds are the BHBs (as determined by color and surface 
gravity) and contaminants to this population are plotted as points- only 1 in 5 points from 
the entire dataset are plotted to avoid overcrowding. The contour lines show the general 
locations, in this color box, of quasars (qso, the reddest in i — z with a locus of about i — z 
= -0.07), main sequence A stars (msa, having a locus at around i — z = -0.12) and white 
dwarfs (wd, the bluest, residing around % — Z—-0.20): other contaminants (oth) are mostly A 
and F stars for which the SSPP failed to assign a gravity or assigned a gravity outside of the 
selection for either BHB or MSA described in §2. In the upper panel, the number and type 
of objects in the BHB selection box (the dashed polygon) are histogrammed as a function 
of g — r color- all of the dataset is presented in this frame. We see that hardly any white 
dwarfs pass this color cut. 
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Fig. 3. — The lower panel shows the distribution of stellar objects (BHBs, MSAs and white 
dwarfs; plotted as dots) in g — r vs g — z color space as compared to that of quasars (plotted 
as circles). Only data which pass the color-color cut shown in Figure 2 are presented here. 
The upper panel shows the number of objects in the BHB selection after this second color- 
color cut as a function of g — r color. Comparison to Figure 2 shows that this second color 
cut is instrumental in removing quasars from the selection sample. 
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Fig. 4. — Color- magnitude diagrams of 10 globular clusters found in the SDSS footprint, 
arranged by radial distance from the sun (Harris (1996) Catalog, 2010 revision values; Table 
1). The drawn boxes indicate the general location (as selected by eye) of the "true" BHBs 
used for the completeness test described in §3. Each box extends 1 magnitude in g from g — r 
= -0.3 to 0.0. Clusters such as NGC 7078 and NGC 5272 show RR Lyrae contamination 
on the red end of the box and clusters such as NGC 6205 and NGC 6341 show sub-dwarf 
contamination on its blue end. 
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Fig. 5. — The northern Galactic hemisphere section of the celestial equator (stripe 10). 
Photometric objects passing the color cuts shown in Figures 2 and 3 are presumed to be 
BHBs and are plotted here. The density of the points are shown via a gray-scaled density 
estimate which is accentuated by the white contours. We see the Hercules-Aquila cloud at 
R.A. w 240° and g P3 16; from R.A. P3 230° to « 200° the Sagittarius dwarf spheroidal's tidal 
stream is apparent at g ~ 19; the Virgo overdensity can be seen at R.A. ^ 200° and q 



7, th e overdensity at R.A. £3 160° and g P3 16 is noted in F-turnoff stars by 



Newberg et al. 



-24 - 



22 




R.A. (deg, J2000) 

Fig. 6. — The southern Galactic hemisphere section of the celestial equator (stripe 82)- 
constructed in a similar fashion to Figure 5. Once again we see evidence for major known 
structures such as Sagittarius at R.A. 30° and g w 18 and the Hercules Aquila cloud at 
R.A « -40° and g « 17. 
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Fig. 7. — The top four frames of this image are SDSS spectra of a typical BHB star, MSA 
star, white dwarf and low redshift quasar, respectively. The botto m frame d e picts the SDSS 



passband functions (ugriz, from left to right) as described by 



Doietal 



feninh and the 



Pan-STARRS passband functions (grizy, from left to right) provided on the Pan-STARRS 



website (http:/ /svn. pan-starrs.ifa.hawaii.edu/trac/ipp/wiki/PSl_Photometric_System). In 
the spectra of the BHB and MSA stars, we see the Paschen features in the 830 nm to 
920 nm range (which falls in the z-band of both systems; expanded to the right)- these 
features will cause a difference in z band based colors for BHB and MSA stars since they 
are susceptible to pressure broadening. 
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Fig. 8. — The left frame shows the Balmer-a line for a composite MSA star (solid line) 
and a composite BHB star (dotted line). Due to Stark pressure broadening, these two 
absorption line profiles have different line widths and shapes which allows separation of 
these two species of star based on their differing surface gravities. The middle frame shows 
the Paschen absorption features in the z band region. These features show a similar behavior 
with the BHB spectra (dotted line) having slimmer absorption lines than the MSA star (solid 
line). Since the continuum at this point is actually a quasi-continuum formed by the Paschen 
features flowing into each other, the BHB star appears to have a higher continuum than the 
MSA star (since broader features would push the effective continuum lower)- this is the 
effect responsible for the differing z band magnitudes. The rightmost frame depicts the 
difference in equivalent width versus core intensity for the Paschen feature residing at 8596A 
(the Paschen 14 line) for BHB stars (crosses; having 3.0 < log(g) < 3.75) and MSA stars 
(circles; having 3.75 < log(g) < 5.0). This plot was constructed using the high signal to noise 
spectra contributing the to combined spectra (described in §5). We see a distinct separation 
between these two species of stars in this parameter space; the pseudo-continuum is higher 
in the BHB stars which causes the feature to have a larger equivalent width and a higher 
core intensity. 
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Table 1: Globular Cluster Information. 



Cluster Name 


R.A. 


Dec 


R (kpc) 


E(B-V) 


Rtidal (') 


g box 


NGC 


4147 


12 10 06.30 


+18 32 33.5 


19.3 


0.02 


6.08 


16.5-17.5 


NGC 


5024 


13 12 55.25 


+ 18 10 05.4 


17.9 


0.02 


18.36 


16.25-17.25 


NGC 


5053 


13 16 27.09 


+ 17 42 00.9 


17.4 


0.01 


11.43 


16-17 


NGC 


5272 


13 42 11.62 


+28 22 38.2 


10.2 


0.01 


28.72 


15-16 


NGC 


5466 


14 05 27.29 


+28 32 04.0 


16.0 


0.00 


15.68 


16-17 


NGC 


5904 


15 18 33.22 


+02 04 51.7 


7.5 


0.03 


23.63 


14.5-15.5 


NGC 


6205 


16 41 41.24 


+36 27 35.5 


7.1 


0.02 


21.01 


14.5-15.5 


NGC 


6341 


17 17 07.39 


+43 08 09.4 


8.3 


0.02 


12.44 


14.6-15.6 


NGC 


7078 


21 29 58.33 


+12 10 01.2 


10.4 


0.10 


27.30 


15-16 


NGC 


7089 


21 33 27.02 


-00 49 23.7 


11.5 


0.06 


12.45 


15.3-16.3 



Note. — All Information from the Harris (1996) Catalog, 2010 revision, excepting the g selection box. 
The g selection box extends from g — r = -0.3 to 0.0 in all cases. 
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Table 2: Globular Cluster Test Statistics. 



Cluster Name 


[Fe/H] 


Total BHBs 


Purity 9ri2 


Completeness^,,.^ 


Purity ugr 


Completeness,^,. 


NGC 


5904 


-1.29 


184 


0.96 


0.38 


0.94 


0.53 


NGC 


5272 


-1.5 


233 


0.94 


0.40 


0.89 


0.54 


NGC 


6205 


-1.53 


251 


0.94 


0.47 


0.84 


0.54 


NGC 


7089 


-1.65 


318 


0.94 


0.37 


0.94 


0.36 


NGC 


4147 


-1.8 


60 


0.96 


0.45 


0.94 


0.48 


NGC 


5466 


-1.98 


115 


1.00 


0.84 


0.96 


0.84 


NGC 


5024 


-2.1 


351 


0.96 


0.70 


0.94 


0.70 


NGC 


5053 


-2.27 


28 


0.92 


0.82 


0.96 


0.79 


NGC 


6341 


-2.31 


154 


0.98 


0.34 


0.96 


0.60 


NGC 


7078 


-2.37 


194 


0.88 


0.57 


0.90 


0.58 


Total 




1888 


0.95 


0.51 


0.92 


0.57 



Note. — This table shows the accuracy and effectiveness of the two tests in extracting BHB stars from 
the chosen globular clusters. Metallicities from the Harris (1996) Catalog, 2010 revision, are included and 
the globular clusters are sorted by this property- one can see a slight trend in both tests to more effectively 
select BHB stars at lower metallicities, especially in the area around -2.0 dex. This is most likely an effect 
of our color cut being formulated by the color distributio n of halo BHB s tars spectroscopically observed by 



the SDSS, which have an average metallicity of -2.0 dex (|Xue et al 



20081) 



